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We study the dynamics of one and two qubits plunged in a q-deformed oscillators environment. Specifically 
we evaluate the decay of quantum coherence and entanglement in time when passing from bosonic to fermionic 
environments. Slowing down of decoherence in the fermionic case is found. The effect only manifests at finite 
temperature. 
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o\ ■ 

Q ■ I. INTRODUCTION 

o : 

Open system dynamics is of uppermost importance in the quantum regime where non classical phenomena turn out to be very 
O . fragile with respect to any noise source. The noise effects are often modeled as the result of an interaction of the system with a 
large number of uncontrollable degrees of freedom, i.e. an environment JTj]. Environments can be assumed as to be composed by 
different kinds of particles, for instance oscillators or spin-i. These objects come, under the mathematical point of view, from 
the realizations of two different algebras (the Heisenberg-Weyl algebra and the Lie algebra su(2)) corresponding to fermionic 
y—i I and bosonic commutation relations. These latter can be seen as two limit cases of more general commutation relations involving 
deformed algebras parameterized by one continuous parameter 12. S 0] • 

Our aim is to analyze the qubit dynamics in an environment of oscillators satisfying suitable q-deformed commutation re- 
lations, such that it permits to continuously interpolate between oscillators and spin-^. Actually, we investigate how quantum 
decoherence phenomena changes in passing from bosonic to fermionic environments. We find a slowing down of decoherence 
£H , in the fermionic case. However, this effect only manifests at finite temperature. 

The paper is organized as follows. In Section|II]we present the model. We then derive the master equation in Section[nJ In 
Section [TV] we study the dynamics of a single qubit and we evaluate its coherence decay. We then study the dynamics of two 
qubits and we evaluate the entanglement decay by distinguishing the case of the two qubits in the same environment (Section 
rvt . from that of the two qubits in separate environments (SectionlVlt. Finally, Section|VIl]is for concluding remarks. 

> ■ 

O ' II- THE MODEL 

o . 

. Let us consider a system (qubit) described by the free Hamiltonian 

H s = na Zl (1) 

with ft the qubit frequency and a, cr\ o z operators satisfying the commutation relations 

[a\a\ = Cz , (2) 
[a,a z ] = 2d, (3) 
[a\a z ] = -2(7+. (4) 

They define the su(2) algebra. Furthermore, we consider an environment composed by an infinite (countable) number of oscil- 
lators whose Hamiltonian reads as Q5|, |6[| 

k 

with ujk the frequency of the fc-th oscillator and Ak, A\, Nk operators satisfying the commutation relations 

[N h ,A k ] = -8 hk A k , (6) 
[N h ,A[] = 5 hk A\. (7) 

They define the Heisenberg-Weyl algebra. We are now going to introduce a deformation of this algebra through the so-called 
"quons" commutation relations [|5|] 

A h A\ - qA\A h = S hk , (8) 
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where q E [—1, 1] is the deformation parameter. It allows us to interpolate between fermions (q = — 1) and bosons (q = 1). 
Intermediate values of q 6 (—1, 1) correspond to the so-called "infinite statistics". 

We assume the system interacting with the environment through the following Hamiltonian 

Hj = J2 X k( A l° + A k<r 1 ), (9) 

k 

where X k denotes the coupling constant of the system with the fc-th environment's oscillator. 



III. MASTER EQUATION 



Quite generally, the master equation for the system density operator p can be derived by using the Born-Markov approximation 
0]] . Hence, it can be formally written as 



drTr E {[Hi(t), [H^t - r),p(t) ® P e}}} ■ 



(10) 



where p E is the initial environment density operator and Tte denotes the trace over environment degrees of freedom. Further- 
more, it is 



(11) 



For the choice of the environment Hamiltonian (0, the dynamical equations are formally identical to the undeformed case. 
The reason is that the interaction Hamiltonian Hi (t) reads as follows 



Hi(t) =J2 Xk (A ae ~ l{Wk ~ n)t + ^*o- t e t(w *- n) *) , 
ft 

by virtue of CD), ®, © and ©. Therefore, from (To), we can write 

roc 

p(t) = - / drTiE Y^\ k \iT k i{p{t)), 
Jo k,l 

where 

F M {p{t)) = F k {t) F t (t - t) p(t) ® PE - F k (t) p{t) ® p E Fi{t - t) 

-F k {t- t) p{t) ® p E Fi (t) + p(t) <g> p E F k (i - t) Fi (i) 

and 

F k {t) = A\ae- i{u>k - n)t + A fe erV (a "^ 0)t . 
We now assume an initial thermal state for the environment at temperature T, 

p E = Z^e- H -' T , 

where 



Z = Tr E {e- H ^ T } 



is the partition function. 

In ( fT3l ) we have nonzero terms of the form 



Tr, 



{p E Al(t)Mt-r)} 



1 



;Tl E { CXp 



2y / i J J 



u k [n k ] q 
2T 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 
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Here we have defined 

r 1 1 - qn 
[»]« = — , 

as the q-deformed number. 

Then, neglecting principal values terms, we obtain from ( TT~8b 



(19) 



/ dTVA fc A/Tr B {p B 4(^(t-T)} e ^ T = / dTY^\l([N{Lo k )] q ) E e-^«^ T 
Jo k,l Jo k 

= ^*l(Wo>k)] q )E6(w k -n), (20) 



where 



[N(u k )] q ) E = i^[n fe ],exp 



Z 

it 



[nk\ q 



(21) 



2T 

Moving to the continuum of frequencies for the environment oscillators, we have 

X 2 k ([N(Lu k )} q ) E 6{u; k - fi) -> J dLoA 2 {u)([N{u;)] q ) E 5(Lu - 0), (22) 

k 

where A 2 (w) accounts for the coupling spectrum as well as for the density of states. As usual, we set A 2 (f2) = 7/2 to be the 
damping rate. Moreover, we get the following distribution Q7|, |8|] 

([N] q ) E ee ([N(Q)] q ) E = eSl/ l_ , (23) 
([N + l] q ) E ee ([N(n) + l] q ) E = en/r _ ■ (24) 

In summary, from (1201 . we have 

f dr^A fc A i Tr £ {p is 4(t)A ; (<-r)} e -^- n ) t +'^-^(*-^ = ^([N] q ) E . (25) 
J k,i 

Other nonzero terms in ( fT3l are 



(26) 



/ dr]T A fc A ; Tr £ {p B A fe (i)4(t - r)} e 'K-n)t-^,-(l)(H = |<[# + 

J k,l 

We finally arrive at the following master equation for the reduced system (the qubit): 

P(t) = ~{W] q )B (™ f p(t) ~ 2a^p(t)a + p(t)a^) -|([iV+ l] q ) E (a^p(t) - 2ap(t)^ + p(t)^a) . 

This equation explicitly shows that the effect of the q-deformation is to change the rates of emission, which is proportional to 
([N] q ) E , and the rate of absorption, proportional to ([N + l] q ) E ( see also [8]). Notice that for T = 0, there are no effects 
coming from the deformation, because for N = we simply have ([0] q ) E = and ([l] q ) E = 1; in other words, the nonlinear 
effects introduced by the q-deformation cannot be observed if the environment transitions only concern the vacuum and the 
states with single excitation. 

IV. ONE QUBIT 

Let us consider the operators appearing in Eq.(l27]) and represent them in matrix form in the computational basis {|0), |1)}, 

JS). 



4 




,i , , , , , 1 

0.5 1 1.5 2 2.5 3 

Yt 



FIG. 1: The plot shows the decay of the coherence bi(t) for a qubit in a quon-environment at temperature T/Q = 1, for different values of 
the deformation parameter. In the inset, it is shown the decay of the population a(t) (solid lines refer again to T/Q = 1, dashed line refers 
T = 0). 



and 



P W- I h(t)-Lb 2 (t) l-o(t) (28) 



where a(t), b\(t) and b 2 (t) are real functions of time to be determined. 

Inserting the above matrices into Eq.(l27b we get the following set of differential equations 



dt 

where for the sake of simplicity we have set 



^-a = -2(A + B)a + 2A, (29) 
d {b l + ib 2 ) = -{A + B){bx + ib 2 ), (30) 



A = (j/2)([N] q ) E , (31) 
B = ( 7 /2)([iV+l] 9 ) B . (32) 



(33) 



The solutions of the differential equations d29l . (l30l read 

-2(A+B)t r , . - 

= A + [a(0)g + A(a(0) + e 2 ^ + ^-l) 

6i(t) = 6 1 (0)e-( A+B )*, (34) 
6a (*) = 6 2 (0)e-( A+B )*. (35) 

Figure [T] shows the decay of the coherence (6i(i)) for a qubit in a quon environment at temperature T/Q, = 1, for different 
values of the deformation parameter. From Eq.s (l3TT l. (l32l . (l34l . it follows that the decay of coherence at q = — 1, and any 
finite temperature, behaves as the decay at T = and any q. In the inset, it is shown the decay of the population a(t) (solid 
lines refer again to T/Q = 1, while dashed line refers T — 0). Thus, the fermionic environment gives rise to the slowest decay 
of coherence and population. The decay of quantum coherence becomes slower and slower when passing from the bosonic to 
fermionic environment. 
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V. TWO QUBITS IN THE SAME ENVIRONMENT 



We now assume the system composed by two identical qubits interacting with the same environment. Then the master equation 
can be written as Eq.d27li simply replacing a with <r 1 + cr 2 , that is 

P(t) = ~(W] q )E ((Tio{p(t) - 2a{p(t)a! + p{t)uxa\ + o 2 a\p{t) 



—2a 2 p(t)(T 2 + p(t)<J2crl + oio\p(t) + <T 2 a\p(t) 
-2a\p(t)cr 2 - 2a\p(t)cr\ + p(t)axu\ + p{t)a 2 a\ 

1{[N + l] q ) E (a\oip{t) - 2a lP {t)a\ + p{t)o\o 1 + a\a 2 p{t) 
-2<j 2 p{t)<j\ + p(t)a\a 2 + a\a 2 p{t) + a\atp{t) 
-2(Jip(t)a\ - 2<j 2 p(t)a\ + p(t)a\a 2 + p(t)(x\a\ 



Then, we proceed in the same way as for the single qubit case. That is, we consider the operators appearing in Eq. 
represent them in matrix form in the computational basis {|00), 1 01) | TO) , |11)}, 
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(36) 



and 



(37) 



and 



Pit) 



/ a(t) b x (t)+ib 2 {t) a(t) + Lc 2 (t) 

h(t)-ib 2 (t) e{t) fi(t) + Lf 2 {t) 

C X (t)-LC 2 (t) fl(t)-lf 2 (t) h(t) 

\ dx{t) - id 2 (t) gi(t)-cg 2 (t) h{t)-ii 2 (t) 



1 



di(t) 

9i(t) 
H(t) 
a(t) - 



VLd 2 (t) 
htsa(t) 
V d 2 (t) 
e(t) - h(t) 



(38) 



where a(t), bi(t), b 2 (t), c\(t), c 2 (t), d\(t), d 2 (t), e(t), fi(t), f 2 (t), gi(t), g 2 (t), h(t), ii(t) and i 2 (t) are real functions of time 
to be determined. In terms of these functions, the master equation is written as a set of coupled differential equations. They are 
reported together with their solutions in Appendix lAl 

In this case, a relevant quantity to study is the entanglement between the two qubits. Specifically we consider the qubits 
initialized in one of the four Bell states 



|0 ± ) = 2-^(100) ± |H)), 
|V±) = 2-V2(|oi> ±|10», 

and then we investigate how entanglement decays. 

We use the concurrence as measure of the degree of entanglement H 

C{ P {t)) - max{0, Ax(t) - \ 2 {t) - A 3 (i) - A 4 (t)} , 

where Ai(i)'s are, in decreasing order, the nonnegative square roots of the moduli of the eigenvalues of p{t)p(t) with 



Pit) 



4) {<r 2 



'(*) (^i 



^2 



(39) 
(40) 



(41) 



(42) 



and p*(t) denotes the complex conjugate of p(t). 

The decay of the concurrence is plotted in Figure|2] The qubits are initialized in the Bell states (1391 . For T/fl = 0, the decay 
of the concurrence is independent from the deformation parameter q. For T/tt > we see the phenomenon of entanglement 
sudden death [ 10]. We notice however that the entanglement death time depends on the value of the deformation parameter q. 
In particular, the slowest decay and the longest lifetime of entanglement is evident for the fermionic case q = —1. The same 
happens when the two-qubit state is initialized in |^ + ) = 2~ 1 / 2 (|01) + 1 10) ) (see inset). The decay of concurrence is slower 
and slower when continuously passing from the bosonic to the fermionic environment. On the contrary, the Bell state is 
invariant under the dynamics of (1361 . thus entanglement in this case is totally preserved. 
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FIG. 2: The plot shows the decay of the concurrence for the two-qubit system in a quon-environment. Both qubits are plunged into the 
same environment. The qubits are initialized in one of the Bell states \<j>±) = 2~ 1/2 (|00) ± |11>); dashed line refer to T/fi = 0, solid 
lines to T/fl = 1 and several values of the deformation parameter q. In the inset the two-qubit state is initialized in the Bell state = 
2" 1/2 (|01) + 1 10) ) . The remaining Bell state \ip-} = 2" 1/2 (|01) - [ 10>) is preserved by the dynamics. 

VL TWO QUBITS IN SEPARATE ENVIRONMENTS 

Here we consider each of the two identical qubit interacting with its own environment. Then the master equation is a straight- 
forward extension of Eq. (l27T l. that is 

p(t) = -^{[N] q ) E (aia\p{t) - 2a\p{t)ai + p{t)<na\ + o 2 a\p{t) - 2a\p{t)a 2 + p(t)a 2 4) 

-7}([N + l] q )s [<rlaxp{t) - 2a x p{t)a\ + p{t)a\ui + a\cT 2 p{t) - 2a 2 p(t)4 + p(t)4a 2 ) ■ (43) 

It can be solved with the same method of d36l l. The corresponding differential equations and their solutions are reported in 
Appendix IE1 

Figure [3] shows the decay of the concurrence in time. The qubits are initialized in the Bell states (|39l l. For T/fl = 0, the 
concurrence decay is independent from the deformation parameter q. Also in this case, for T/fl > 0, we see the phenomenon of 
entanglement sudden death fioll . We notice that the entanglement death time depends on the value of the deformation parameter 
q. In particular, the slowest decay and the longest lifetime of entanglement is for the fermionic case q = — 1. The same happens 
when the two-qubit state is initialized in \i()±) — 2 _1 / 2 (|01) ± 1 10) ) (see inset). The decay of entanglement becomes slower 
and slower when passing from the bosonic to the fermionic environment. In this case there is no maximally entangled state that 
remains invariant under the dynamics. 

VII. CONCLUSION 

In conclusion, we have analyzed the qubit dynamics in an environment of oscillators satisfying suitable g-deformed com- 
mutation relations, such that it permits to interpolate between oscillators and spin-^ particles. Specifically we have evaluated 
the decay of quantum coherence and entanglement in time when passing from bosonic to fermionic environments. The general 
behavior is that, at finite temperature, coherence and entanglement decay slower and slower when continuously passing from 
bosonic to fermionic environments. 

Our work sheds further light on the mechanism of loosing quantum coherence and paves the way for a deeper algebraic 
analysis of this phenomenon. Moreover it could be useful for describing realistic physical situations where the assumption of an 
interaction with an environment of solely oscillators (resp. spin-i) particle turns out to be oversimplified. 
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FIG. 3: The plot shows the decay of the concurrence of the two-qubit system in a quon-environment. Each qubit is subject to independent and 
identical environment. The qubits are initialized in one of the Bell states \<f>±) = 2~ 1/2 (|00) ± |11»; dashed line refer to T/Q = 0, solid 
lines to T/fi = 1 and several values of the deformation parameter q. In the inset the two-qubit state is initialized in one of the Bell states 

|V.±)=2- 1 / 2 (|01)±|10)). 



APPENDIX A: TWO QUBITS IN THE SAME ENVIRONMENT 

Using the parametrization in d38l ). the master equation d36*l i translates in the following set of differential equations 

^-a = -ABa + 2A(e + 2/i + h) (Al) 
— (h + ibi) = -B(3b! + 3ib 2 + ci + lc 2 ) 



dt 

d 
di 

d 



-A{bi + ib 2 + ci + LC2 - 2. 9l - 2tg 2 - 2z x - 2« 2 ) (A2) 
(ci + lc 2 ) = —Bipx + ib 2 + 3ci + 3ic 2 ) 

-A(6i + ib 2 + a + lc 2 - 2.9i - 2ig 2 - 2h - 2w 2 ) (A3) 

Adt + uh) = -2(A + B)(d 1 +id 2 ) (A4) 
at 

-f-e = 2B{a - e - /i) + A(a + 2e + /i + /i - 1) (A5) 
at 

+ - B{2a-e-h-2f 1 -2tf 2 ) 

dt 

d_ 

di 

d 



+A(2 - 2a - 3e - 3/i - 2/i - 2t/ 2 ) (A6) 
(.9i + ^2) = S(26i + 2ife 2 + 2ci + 2tc 2 - g\ - tg 2 —ii — m 2 ) 

-^(3^1 + 3tg 2 + ii + m 2 ) (A7) 

h = 2B(a -h-fi)- 2A(a + e + 2h + /1 - 1) (A8) 



— (ii + w 2 ) = B(2bi + 2ib 2 + 2ci + 2ic 2 - 51 — tc/ 2 — i% - d 2 ) 
dt 

-A( 9l + ig 2 + Mi + 3w 2 ) (A9) 
The solutions of the differential equations (IAlh -( IA9b read 
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hit) 



-4(A+s/AB+B)t 



{-10a(0)AB 2 e iM (e 8 ^ - l) 

) (e(0) + 2/i(0) + fc(0)) 



8VAB(A 2 -3AB + B 2 ) 

+2a(0)B 3 e 4At ( e 8 ^4i?t _ ^ + A ^ e ^t ^sSTB, 
+A 2 Be 4At ( e 8 ^4B* - l) (4a(0) - 3e(0) - 6/1 (0) - 3/i(0)) 
-2VI^B [ 2e 4 ( A+v ^ +B ) t - a(0) (e 4At - 2e< VAB+B ) t + e ^t+sVABt^ 
+2e 4 ( yZg+B ) t (e(0) + ft(0) - 1) 



e 4At (e(0) + 2/ 1 (0) + /i(0))-e 



(e(0) + 2/i(0) + /i(0)) 



+VAB5 r_ 4e 4(A+VAB+B)t _ 2e 4 (^4iJ+ B ) t ( e (0) - 2/i(0) + /i(0) - 2) 
+e 4At (e(0) + 2,A(0) + h(0)) + e 4At+sVABt (e(0) + 2/i(0) + h(0)j 

+ VA^ [l 2e < A +^B+B)t _ 6a(0 ) ^At _ 2e i(^AB+B)t + e AAt+ S VABt^ 

-5e 4At (e{0) + 2A(0) + h(0)) - 5e 4A *+ 8v ^*(e(0) + 2/i(0) + h(0)) 
+2e< VAB+B ) t (5e(0) - 2/i(0) + 5/i(0) - 6)] } , 



(A15) 



/ 2 (i) - / 2 (0)e- 2 (' 4 + B ) t , 



(A16) 



~(3A+3B+0)t 

= 492 {i?[4e(6 1 (0) + Cl (0))(-l + e 2et ) + 14A( ffl (0) + n(0)) 

-9 ( 5l (0) + i!(0)) + 6e 2et ( 5l (0) + ix(0)) 
+28^ +3S + e ) t ( 5l (0) - ix(0)) + 14Ae 2et ( 5l (0) + h(0))" 

+B 2 [(l + e 2et + 2e {A+3B+e ^) 5l (0) + (l + e 2et - 2^ A+3B +^ t ) h(0) 

9 (1 - e 2et ) ( 5l (0) + ix(0)) + A ((l + e 2et + 2e ^+ 3B + e )>) ffl (0) 



+i4 



+ (l + e^-2e< A+3B+e >*) *!(()))]} 



(A17) 
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h(t) 



-2(3A+2v / AB+3B)t 



[-2a{0)AB 2 e 2{3A+B ^ (e 



lABt 



8VAB {A 2 - 3AB + B 2 ) 
+2a(0)B 3 e 2 ^ A+B '> t ( e 8 ^ _ i) 

_ A 3 e 2( 3 A+B)t ^sVABt _ {j (e(Q) + 2/i(0) + m) 

+3A 2 Be 2( - 3A+B '> t ( e sVABt - l) (e(0) + 2^(0) + h(0)) 
+VAB^e 2 ^ A+B ^ U e <A+VAS+B)t _ 4e 2(VA+VB) 2 t^ _ h(Q)) 



+2e 



4VABt+ABt 



(e(0)-2/ 1 (0) + / l (0) 



+e AAt (e(0) + 2A(0) + h(0)) + e AAt+8VABt (e(0) + 2/i(0) + /i(0)) 
-VA 3 B 3 e 2(A+B)t \i2 e A ( A +^ B+B ) t 



e 4At _|_ e 4At+sVABt 



-6a(0) (■ 

-12e 2 ( v ^ +v ^) 2 *(e(0) - /i(0)) + e 4At (e(0) + 2/i(0) + /i(0)) 



4\/ABt+4sA 



+e 



4At+8VABt 



(e(0) + 2/i(0) + /i(0)) 



+2e 



A^ABt+ABt 



(5e(0)-2/i(0) + 5ft(0)-6) 



+ 2VI^B [-a(O) (e 2(iM+B)t + e 2(3A+W3B+u)* _ 2e 2(A+2VAB+3S)t^ 
+ 2 e 2 ( A + 2 ^B+ 2B ) t ^ e 2 (2-4+-B)t 

+ e 2S *( e (o) + Mo)-i))]}, 



e 2At (e(0)-MO)) 



(A18) 



ii(t) 



e 



3 -(3A+3B+e)t 



49 2 



{B [49 (6i(0) + ci(0)) (-1 + e 



2(-)n 



+14A ( ff i(0) + u(0)) - 9 (0) + *i(0)) + 6e 2e * ( Sl (0) + u(0)) 
-28Ae^ A+3B+@ ^ ( 5l (0) - *i(0)) + 14Ae 2et ( 5l (0) + ii(0))" 

2 + e 26* _ 2e (A+3B + e )t ^ gi(0) + ^ + e 20* + 2e (A+3B + e)t^ 



+B 
+A 2 



(l + e 2e * - 2e ( A + 3B+e ) t ) 3l (0) + (l + e 2e * + 2e ( - A + 3B+e ^ h(0) 



+A [9(l- e 2et )( 3l (0) + Jl (0))]} 



(A19) 



where 9 = V 'A 2 + UAB + B 2 . The other solutions b 2 (t), c 2 (t), d 2 (t), g 2 {t) and i 2 (t) can be obtained from bi(t), a(t), 
d\(t), gi(t) and i\{t) respectively by simply replacing the subscripts 1 — > 2. 



APPENDIX B: TWO QUBITS IN SEPARATE ENVIRONMENTS 



Using the parametrization in 



the master equation (l43l l translates in the following set of differential equations 

d 



dt 



dt 
d 

~dt 

d 



dt 
d 

dt 



dt 



-4Ba + 2A(e + h), 
(61 + bb 2 ) = -3B(h + ib 2 ) - A(h + ib 2 - 2ii - 2w 2 ), 
-W(ci + lc 2 ) - A(c\ + lc 2 ~ 2gx - 2ug 2 ) 

-2(A + B)(di + ida), 

2B(a - e) - 2A(a + 2e + h- 1), 

(ji + ih) = -2(A + B)(h+if 2 ), 

(gi + ig 2 ) = B(2c\ + 2lc 2 - g-y - ug 2 ) - 3A(#i + Lg 2 ), 
d 

2B(a -h)- 2A(a + e + 2h- 1), 
B(2bi + 2ib 2 - h - ii 2 ) - 3A(ii + ii 2 ). 



(ci + LC 2 ) 

(di + td 2 ) 
d 

— e 
dt 



dt 



[h + "2) 



The solutions of the differential equations (IB 1 b - (IB9b read 

e -A(A+B)t 



a(t) = 



(A + Bf 



|o(0)S 



1 , ^2 
2(A+B)t _ 



Cl (t) 
di(t) 
e(t) 



e -3(A+B)t 

A + 5 

e -3(A+B)t 

A + 5 

e -i(A+B)t 



3(0) (2^+^ - l) 

l) (e a ( A+B )* + e(0)+/i(0)-l)" 
2a(0)e 2 ( A+s > i + (e(0) + h(0)) ( e 2 ^+ s )' - l)] Ai?} , 

{6i(0)fl + A [6 1 (0)e 2 <' 4 + i3 ) t + (e 2 ( A+J ^ - l) ^(0)1 } , 

{ci(0)B + A [c 1 (0)e 2 ( A + B ) t + (e 2 ^+ s )* - l) 5l (0)] } , 



(A + Bf 



{ [a(0) (l - e 2 ^^) + e(0) - e 2 ( A + s )* (ft(o) - 1) + h(0) - 1 



e(0) + e 4 ( A+s )* + e 2 ^ 5 )* ( e (0) - fc(0) - 1) + ft(0) 
a(0) (e 2 ^+ B )* - l) + e(0)e 2 ( A+B >*] £ 2 } , 



AS 



fc(t) 



h(Q)e-^ A + B ^ 

e -3(A+B)t 



A + B 

e -i(A+B)t 

(A + B) 2 



{ 91 (0)B + A [ 9l (0)e 2{A+B)t + (e« A+B » - l) c x (0)] } , 
{ [o(0) (l - e 2 ^)*) - e 2(A+B ^ (e(0) - 1) + e(0) + h(0) - 1 



A 



e(0) + e 4 ( A + B )* + /i(0) 



,2(A+B)t 



(ft(0) - e(0) - 1) 



AB 



o(0) ( 



= 2(A+B)t 



-3(A+B)t 

A + B 



{t^B + A ^(O)^^ + (e 2{ - A+B ^ - l) 6 X (0)] } 
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The other solutions 62(^)1 C2(t), dsfy), f2(t), g2(t) and t2{t) can be obtained from &i(t), ci(t), 51 (t) and 

respectively by simply replacing the subscripts 1 — > 2. 



[1] C. W. Gardiner, Quantum Noise, Springer, Berlin 1991. 

[2] V. Kuryshkin, Annales de la Fondation Louis de-Broglie 5 (1980), 111. 

[3] P. Kulish and E. Damaskinsky, J. Phys. A 23 (1990), L415. 

[4] A. J. Macfarlane, J. Phys. A 22 (1989), 4581; 

L. C. Biedenharn, J. Phys. A 22 (1989), L873; 

C.-P. Sun and H. C. Fu, J. Phys. A 22 (1989), L983. 
[5] O. W. Greenberg, Phys. Rev. Lett. 64 (1990), 705; Phys. Rev. D 43 (1991), 41 1 1. 
[6] R. N. Mohapatra, Phys. Lett. B 242 (1990), 407. 
[7] M. Chaichian et al., J. Phys. A 26 (1993) , 4017. 
[8] J. W. Goodison and D. J. Toms |arXiv:hep-th/9410096l 
[9] W. K. Wootters, Phys. Rev. Lett. 80 (1998), 2245. 
[10] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93 (2004), 140404; 

A. Al-Qasimi and D. F. V. James, Phys. Rev. A 77 (2008), 012117. 



